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We have measured the distributions of protons and deuterons produced in the 20% most central 
Au+Au collisions at RHIC (a/s^ — 200 GeV) over a very wide range of transverse and longi- 
tudinal momentum. Near midrapidity we have also measured the distribution of antiprotons and 
antideuterons. We present our results in the context of coalescence models. In particular we ex- 
tract the "homogeneity volume" and the average phase-space density for protons and antiprotons. 
Near central rapidity the coalescence parameter B2{pt) and the space-averaged phase-space density 
{f){pT) are very similar for both protons and antiprotons. For protons we see little variation of 
either B2{pt) or the space-averaged phase-space density as the rapidity increases from to 3. How- 
ever, these quantities depend strongly on pt at all rapidities. These results are in contrast to data 
from lower energy collisions where the proton and antiproton phase-space densities are different at 
j/=0 and both B2 and (/) depend strongly on rapidity. 



PACS numbers: 25.75.Gz, 25.75.Ld, 13.85.2t, 25.40.Ve 

I. INTRODUCTION 

Deuterons detected in heavy ion collisions are conven- 
tionally thought to be produced predominantly via a pro- 
cess called coalescence. Protons and neutrons that are 
close enough in phase-space, (i.e. in position and mo- 
mentum space) "coalesce" to form deuterons. The sur- 
rounding medium created in the A-l-A collisions enables 
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this "2 ->1" process to proceed while conserving energy 
and momentum. There is a wide range of evidence to sug- 
gest that a dense system of strongly interacting partons 
is created in heavy ion collisions at Relativistic Heavy Ion 
Collider (RHIC) energies Hj. In the hot and dense system 
produced in high-energy ion collisions, the coalescence of 
nucleons into deuterons cannot even begin before the the 
partons have frozen out into hadrons. Even then, the low 
binding energy of deuterons (^2.2 MeV), ensures that 
only deuterons produced late will survive to reach the 
detectors. Thus our sample is dominated by deuterons 
produced close to thermal freeze-out, where the nuclear 
density is low [H [3] . 

Early coalescence models assumed that the phase- 
space density of clusters is proportional to the product of 
the phase-space densities of individual nucleons that coa- 
lesced into them and that the momentum of the deuteron 
is trivially the sum of the miclcon momenta [5^45] . More 
recent coalescence models add insight into the nature of 
the phenomenon, but still relate cluster production to 
the product of phase-space densities [7j- For deuterons 
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where p is the momentum of the proton and B2 is de- 
fined as the coefficient hnking the square of the nucleon 
distribution to that of the proton distribution. 

Ideally for a coalescence study one would measure the 
spectra of both protons and neutrons. At ^/s^ = 4.9 
GeV the n/p ratio measured in the 10% most central Au 
on Pb collisions has a value of 1.19 ± .08 and is indepen- 
dent of rriT in a wide range of rapidity (1.6 to 2.4) [5], 
whereas the same ratio for the incident nuclei is equal 
to 1.52. (Measurements by other experiments imply that 
most of the isospin lost by the neutrons goes into an ex- 
cess of 7r~ over 7r+ [3[in]) At higher energies we expect 
n/p to be close to 1 over a wide rapidity range. 

In thermodynamic models that assume thermalized 
distributions of nucleons, B2 carries information about 
the effective volume (in coordinate space) of the nucle- 
ons when they coalescence. This is also known as the 
homogeneity volume, and is defined as that volume over 
which nucleons, are close enough in momentum that they 
can coalesce [IHS]. At RHIC the homogeneity volume of 
the emitting source has been studied in great detail at 
central rapidity but few measurements are available in 
the forward region. Consequently the deuteron coales- 
cence analysis gives BRAHMS a unique opportunity to 
study the volume of the proton source over a wide range 
of rapidity. 

The paper is organized as follows. This introduction 
continues with a short review of the homogeneity volume 
and the phase-space density averaged over coordinate 
space. Sec. II details the detector setup, the selection 
of the central events used for this analysis, the tracking 
algorithm, particle identification and the feed down cor- 
rections applied to the proton and antiproton spectra. In 
Sec. Ill we describe the fully corrected spectra and the 
inferred values of B2 and average phase-space density as 
functions of transverse momentum and rapidity. Finally, 
we summarize our results and the physics that we extract 
from them. 
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where rrip denotes the proton mass. The left-hand side 
of this equation just represents the convolution of the 
proton source size and the Gaussian representation of the 
deuteron wave-function. For small systems 5 dominates 
over Rq and the measurements have only little sensitivity 
to the proton source size. One advantage of this ansatz 
is that it facilitates comparison to interferometry radii. 
We show in Sec. Ill that this model works well for our 
data since we find consistency between our estimate of 
the proton source size and Hanbury Brown- Twiss (HBT) 
results. At lower beam energies Rq has been found to be 
very consistent with radii measured from tttt, KK and 
pp correlations [T3] • The deuteron wave- function is more 
accurately represented by the Hulthen form, which has 
an exponential tail [14]. For the B2 values reported in 
this paper the error in Rq from using Eq. [2] is less than 
0.2 fm [13]. 

For a variety of ion beams ranging from He to Ar 
incident on nuclear targets at BEVALAC energies |15) . 
as well as p+K fixed target experiments at FNAL [16j . 
KEK [T7] and SPS [T^, the measured B2 values are 
independent of energy and pT and are consistent with 
measurements of the spatial parameters of the deuteron 
wave-function. At higher A-f A energies including AGS 
Au-hAu fixed target ,Js^=4,.'d GeV [MH], SPS S-^S, 
S-HPb and Pb-HPb at ys^« 18 GeV [MS], and 
Au+Au {,/s^=2m GeV) at the RHIC colhder [MSQ], 
B2 decreases with energy and increases with pT- This is 
consistent with the formation of deuterons in an expand- 
ing medium [31H33| . Near midrapidity B2 decreases by 
a factor of « 20 as the center of mass energy increases 
from y/s^Ki 1 to 17.3 GeV [26] before flattening out at 
RHIC energies [M', 'HS] . This is similar to the behavior 
seen in HBT radii [36i l37j . 



B. Average Phase-space Density 



A. B2 and the Homogeneity Volume 

The link between B2 and the homogeneity volume can 
be seen in coalescence models that use the so-called sud- 
den approximation technique to write the phase-space 
density of deuterons as an overlap of the nucleon wave 
functions with the deuteron wave function. The volume 
of the nucleon system is introduced through the normal- 
ization of their wave functions jllj . One can make a 
reasonable approximation of the deuteron's wave func- 
tion by assuming it is a Gaussian of width (5 = 2.8 fm. 
If we further assume that the region where coalescence 
occurs also has a Gaussian spatial profile with width Rq^ 



The average of the phase-space density /(x, p) over the 
system volume at freeze-out time is a quantity, which, 
when compared to Bose-Einstein or Fermi-Dirac statis- 
tics, allows the assessment of the degree of chemical or 
kinetic equilibrium reached by the system at that stage 
[381 139] . This quantity also carries information about 
the possible multiple-particle symmetrization effects i.e. 
pion condensates. Values greater than 1.0 would indi- 
cate the presence of quantum effects |38[ . This quantity 
is a count of the states that the system can occupy, and 
as such, is used to define the contribution of particular 
types of particles to the overall entropy of the system. 
The spatial average of the phase-space density is defined 
as the following ratio: 
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where the integration is carried over spatial coordinates 
bound by the homogeneity volume of the system at 
freeze-out. The formal definition of phase-space density 
/(x, p) for a particle of spin J is written as: 
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For a system in chemical equilibrium at a temperature 
T and chemical potential fj,, 
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where E is the energy and ±1 selects bosons or fermions 
respectively. For a dilute system, (i.e. / ^ 1), Eq. [s] 
gives: 
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Since Ed — mTCOsh{y), one would expect the phase- 
space density to be an exponential in rriT- Note in this 
simple derivation, we are ignoring the collective motion of 
the of particles. At y/s^ = 17.3 GeV, it was found that 
strong longitudinal flow could significantly reduce the 
space-averaged phase-space density of pions [10]. Also 
at this energy the inverse slope of the phase-space den- 
sity was found to increase with particle mass in a man- 
ner suggestive of transverse flow [IT] . Using the fact that 
the deuteron has energy Ed = En 
P = 2p, Eq. [6] implies that 



Ep and momentum 



/rf(x,P) = /p(x,p) • /„(x,p) = /p2(x,p). (7) 

To extract the average phase-space density of protons 
one can then replace the square term in the numerator 
of Eq. |3]by the phase-space density of deuterons: 
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Furthermore, one can make use of the assumption that 
deuterons are formed by coalescence and satisfy Eqs. [l] 
and[2Jto obtain an expression for the average phase-space 
of protons similar to what Bertsch originally suggested 
for pions. Dividing the spectrum by the product of the 
HBTradn gives [38l[39]: 
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The phase-space densities calculated using Eq. [8] have 
the expected exponential dependence in ttit. When using 



Eq. [9] one introduces the assumption that deuterons are 
produced via coalescence and that the homogeneity vol- 
ume extracted from deuteron distributions is also the vol- 
ume of the proton source. The calculation of the space- 
averaged phase-space density of protons using these as- 
sumptions is close in value to the one obtained from an 
assumed thermal equilibrium and has the same behavior 
in rriT- This suggests that the system has evolved into 
chemical equilibrium at freeze-out. While chemical equi- 
librium has been very well established at RHIC energies 
for the system at midrapidity, these data allow us to test 
if the matter at forward rapidity also reaches equilibrium. 
To ease the comparison with lower energy data we have 
decided to use Eq. |8] to calculate the phase density. 

The coalescence parameter B2 and the space-averaged 
phase-space density recast the information contained in 
the proton and deuteron spectra into "dynamic" and 
"chemical" terms. The coalescence parameter B2 can be 
interpreted in terms of the homogeneity volume, which 
depends upon the temperature of the system and the ra- 
dial flow. Indeed one confirmation that we are actually 
seeing coalescence would be to see if this radius was con- 
sistent with the appropriate HBT radii. This was checked 
at y^s^=17.3 GeV by comparing Rq to radii extracted 
from TTTT, KK, and pp correlation measurements [1311421- 
I44| . In this paper, we use the very large angular and 
momentum ranges of the two BRAHMS spectrometers 
to measure the rapidity dependence of the homogeneity 
volume and space-averaged phase-space density of the 
(anti-)proton distributions for central Au-|-Au collisions. 



II. ANALYSIS 

The data in this paper were collected by the BRAHMS 
experiment during 2004. We present proton and deuteron 
spectra at ^/s^ =200 GeV AuAu collisions with a cen- 
trality range of 0-20%. The statistics for deuterons at 
high rapidity for more peripheral data is limited, so a cen- 
trality dependence analysis was not warranted. The data 
are analyzed in four rapidity bins: (—0.1,0.1), (0.5, 1.2), 
(1.5, 2.5), and (2.8, 3.2). We have also measured antipro- 
ton and antideuteron spectra at ?; sa and y w 0.8. 



A. Detector System 

The BRAHMS experimental setup consists of two mov- 
able magnetic spectrometers, the Forward Spectrome- 
ter (FS) that can be rotated from 2.3° to 15°, and the 
mid-rapidity Spectrometer (MRS) that can be rotated 
from 90° to 30° degrees relative to the beam line. Fast 
Cherenkov counters at high rapidities were used to mea- 
sure luminosity, to determine the interaction vertex, and 
to provide a start time for time-of-flight measurements. 

The MRS has two time projection chambers (TPCs), 
TPMl and TPM2, situated in field free regions in front 
and behind a dipole magnet. This assembly is followed 



by a highly segmented scintillator time-of-flight (TOF) 
wall at 4.51 m. The FS consists of 4 dipole magnets 
Dl, D2, D3 and D4 with a total bending power of up 
to 9.2 Tm. The spectrometer has 5 tracking stations Tl 
through T5. Tl and T2 are TPCs placed in front of 
and after the second dipole D2. T3, T4 and T5 are drift 
chambers with T3 in front of D3, T4 between D3 and 
D4, and T5 after D4 and just in front of the ring imag- 
ing Cherenkov (RICH) [^. Details on the BRAHMS 
experimental setup can be found in Ref. |46j . 



B. Centrality selection 

The centrality of the collisions at the BRAHMS inter- 
action region was extracted from charged particle mul- 
tiplicity densities measured with a multiplicity array 
(MA). The MA array consists of a coaxial arrangement 
of Si strip detectors and scintillator tiles surrounding the 
intersection region [37] . The pseudo-rapidity coverage of 
the MA is approximately —2.2 < ?/ < 2.2. The centrality 
selection is derived from minimum-bias trigger events, 
which are defined using two zero degree calorimeters 
(ZDC), requiring energy deposit equivalent to at least 
one neutron in each of the two detectors and also re- 
quiring a signal in the MA in order to reject Coulomb 
dissociation events [JSl HH] . The uncertainty in the cen- 
trality determination from our multiplicity distributions 
is estimated to be ±4% of the size of the bin for the 0%- 
20% bin. The fraction of the inclusive yield lost by the 
minimum-bias trigger is estimated to be about 4% and 
was corrected for. 



C. Tracking 

Track reconstruction starts by finding local track seg- 
ments in the TPCs and (for the forward spectrometer) 
in the drift chambers. These detectors are in field free 
regions and so these local tracks are straight lines. Lo- 
cal tracks on either side of a magnet are matched using 
the effective edge approximation generating what we call 
locally matched tracks. The locally matched tracks are 
then combined in the FS to form complete tracks. The 
complete tracks are refitted to extract the best measure- 
ment of the particle three momentum. Tracks in the FS 
are required to project through the magnet Dl onto the 
nominal beamline. Track quality cuts are applied for the 
final track selection. The momentum resolution at full 
field ( where all data presented in this paper were taken) 
is: 5p/p = O.OOOSp. 

The momentum distributions are corrected for tracking 
detector efficiency and geometrical acceptance. Tracking 
and matching efficiencies for each of the 5 tracking sta- 
tions in the FS were calculated by constructing full tracks 
using just four track segments and evaluating the effi- 
ciency in the fifth station by comparing the predicted po- 
sition and direction of the interpolated or extra-polated 



full track in that station with the known local segments. 
The local track efficiency was evaluated as function of 
spectrometer angle and field setting, as well as position 
and direction of the local track segments. The overall 
tracking efficiency is typically 60%, and is known to a 
few percent. In the MRS the tracking efficiencies are de- 
termined from Monte Carlo simulation using embedding 
techniques. The systematic error is below 4%. We cor- 
rect for absorption, multiple scattering, and energy loss 
in the detectors using GEANT 3.21 [49]. The magnitude 
of these corrections on the particle yields depends on the 
particle momenta and the spectrometer positions, and in 
the FS it is about 20% with an estimated uncertainty of 
2%. In the MRS the corrections are less than 8% with 
a small systematic error. We assume that deuteron ab- 
sorption factorizes from other effects and model it as the 
square of the proton absorption correction (at pt/2) as 
done in [T3l. More details about tracking in the MRS 
can be found in Ref. [501. 



D. Particle Identification 



In the midrapidity region, particle identification is 
done using time of flight measurements, while at forward 
rapidities we use the RICH [J^. At y=0 and 0.8 we 
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select (anti)protons using a 3cr cut on m — p ( 
where L is the length of the path followed by a particle, t 
its time of flight, and c the speed of light. For dcuterons 
a simple 3cr cut around the (anti) deuteron mean mass is 
used. These deuteron yields are corrected for background 
that arises from false matching of tracks with background 
hits in the TOF wall. This background is determined 
from data in the mass region outside the deuteron peak, 
and contributes to the systematic error of the deuteron 
yields by ^ 5 — 10% decreasing with pT- At forward 
rapidities the RICH provides proton/deutcron FID sepa- 
ration. The ring radius (R) in the RICH depends on the 
particle's velocity. For a particle of mass to in the RICH 
R w 9.2(1- (p/pthres)'^) cm where pthres = 16.1-mc. The 
FID performance is shown in Fig.[l] where the bands rep- 
resent the limits used for the analysis. Because of inef- 
ficiencies about 3% of the particles moving faster than 
their Cherenkov threshold don't produce an identified 
ring. Protons and deuterons with momenta between the 
kaon and proton thresholds do not make a ring. The 
contamination from mis-identified kaons and pions is sub- 
tracted to deduce the proton spectrum in this momen- 
tum range (neglecting the small fraction of deuterons). 
The contamination correction to the proton spectrum is 
-- 10% at pt=0.5 GeV/c and drops to ~ 2% at pt=1-5 
GeV/c. The error on this correction is neglible. This pro- 
cedure cannot be used in the region between the proton 
and deuteron thresholds since the contamination from 
pions, kaons, and protons dominates the small deuteron 
yield. 
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FIG. 1. (color online) Ring radius versus momentum for par- 
ticles at 3/ ~ 3 showing the pion, kaon, proton and deuteron 
separation. The bands show the proton and deuteron selec- 
tion used in the analysis. 
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E. Feed-down corrections 

We have corrected our data to account for the hyperons 
that decay into protons using the method described in 
[5^ . The correction factor, C, is given by: 






(10) 



where Np is the number of primary protons and A'a and 
and N-^+ are the number of protons coming from A and 
S+ decay, respectively. Near central rapidity we have 
used the A spectra measured by PHENIX [51] and es- 
timated the S/A ratio from lower energy measurements 
|52j . Since there are no measurements of As at forward 
rapidities we have estimated the A/p ratio based on ther- 
mal models that were fitted to the rapidity densities of 
charged pions, kaons, protons and antiprotons measured 
by BRAHMS in the forward region [53]. 

With our model calculations, we find that the ratio 
^j^/-g-^ varies slowly with rapidity up to rapidities 
y sa 4. The systematic error from uncertainties on the 
yields and the model extrapolation is estimated to be 
less than 3%. The correction factor also depends on the 
Pt dependence of the A/p ratio. BRAHMS has found 
that the mean transverse kinetic energy scales linearly 
with the mass of the particle with a slope that depends 
only weakly on rapidity [51j. We have used these sys- 
tematics to estimate the inverse slope (T) of the A ttit- 
distribution. To estimate the systematic error on the px 
dependence of the correction factor we have taken the 
limiting cases of T\ — Tp and T\ — Tp ■ rriA/rnp. This 
produces an error on the correction factor that is almost 
zero at Pt — I GeV/c and reaches —9%, -1-6% at pr 
= 2 GeV/c. The correction factor as a function of pt 
and rapidity is shown in Fig. [2] It varies only weakly 
with rapidity and has a systenratic error that is small in 



FIG. 2. (Color online) pr dependent feed-down correction 
factors for protons versus pr and rapidity. The three sets 
of lines show the effect of various scenarios for the difference 
in shape between proton and A spectra. Systematic errors 
on the correction due to uncertainties in the yield ratio of 
protons and lambdas are shown by the bands to the left of 
the plots. The curves for antiprotons are very similar. 



comparison to the statistical errors on B2 and the space- 
averaged phase-space density. At p^ = 2 GeV/c the to- 
tal error from the feed-down correction reaches a value of 
19% for B2 and 10% for the space-averaged phase-space 
density. 



III. RESULTS 

Figure [3] shows the invariant proton and deuteron spec- 
tra versus px for the four measured rapidity bins. The 
systematic error on the normalization of the spectra is 
estimated to be about 5% with a point-to-point error in 
the spectra of ~ 8% arising from the merging of spec- 
trometer settings and FID and uncertainties in tracking 
efficiency determination. We fit exponentials in rny to 
the spectra and extract the invariant yields, dN/dy and 
inverse slopes T. These are listed in Table |lj About 40- 
50% (dependent on rapidity) of the yield of the deuterons 
is contained within our acceptance. For protons the frac- 
tions are slightly higher 55-60%. The proton yield falls 
slightly fronr y=0 to y=3, but no firm conclusion can be 
made for the deuterons due to the limited px range at 
high rapidity. The inverse slopes tend to decrease with 
rapidity, suggesting a decrease of radial flow. The inverse 
slopes of the deuterons are somewhat higher than those 
of the protons. The average over all rapidities for the 
ratio of the inverse slopes is 1.6 ± 0.1. For antideuterons 
the inverse slopes are a factor of 2.0 ± 0.2 higher than 



y 


PT fit 


Proton 

dN/dy 


TMeV 


PT fit 


Deuteron 

dN/dy 


TMeV 


0.0 


0.7-4.0 


27.9 ±0.1 
20.8 ±0.1 


354 ±2 
352 ±1 


1.5-3.3 


0.093 ±0.008 
0.033 ±0.004 


570 ± 70 
870 ± 220 


0.8 


0.5-4.0 


26.0 ±0.1 
17.9±0.1 


356 ±1 
361 ±1 


1.5-3.3 


0.068 ±0.003 
0.031 ±0.002 


610 ±30 
700 ± 80 


2.0 


0.5-3.0 


20.9 ±0.2 


314 ±3 


0.8-2.8 


0.082 ±0.011 


460 ±110 


3.0 


0.5-4.0 


23.4±0.1 


282 ±1 


2.5-3.5 




310 ±90 



TABLE I. Proton and deuteron yields, dN/dy, and inverse 
slopes, T (MeV), derived from fitting spectra. The lower rows 
at rapidities y = Q and y = 0.8 are for the anti-particles. The 
errors listed are statistical only. For dN/dy the systematic 
errors due to extrapolation to low pt is 5-10% with an addi- 
tional 8% from normalization and various corrections. The 
systematic errors on the inverse slope are less than 10%. 



that of antiprotons. 

Figure |4] shows B2 versus pt and rapidity. B2 increases 
with Pt, which is consistent with previous experiments 
[ini 135] ■ Using Eq. |2J we find that at central rapidity the 
source radius Rg falls from 4.2 ±0.2 fm to 3.1 ±0.4 fm as 
rriT increases from 1.2 to 1.9 GeV/c^. This is consistent 
with the rriT dependence of HBT radii, R ex 1/ ^/rnr that 
has been observed by PHENIX ^ and STAR [^ and 
also at SPS energies [43]. The sohd hne in each rapidity 
panel represents an exponential fit to our data at y=Q. 
We see no evidence of any rapidity dependence of B2{pt)- 
The proton and antiproton B2 values are very close at 
this energy implying a similar source size. This is in 
contrast to the measurements at ^s^ = 17.3 GeV where 
the antiproton source volume was found to be somewhat 
larger than the proton source volume [T5] . 

Figure [5] shows the average phase-space density 
(/)(wt) for protons and antiprotons as a function of ra- 
pidity. These values are calculated using Eq. |8] The 
space-averaged phase-space density decreases as the uit 
increases as expected from Eq. [5] The solid curve in each 
panel of Fig. [s] is an exponential fit to the proton density 
at y=0. We see little rapidity dependence of {f){pT)- 

From Eq. [5] we would expect the ratio of the proton 
and antiproton phase densities to be flat as a function 
oi Pt- Fitting a constant to {fp)/{fp) yields a x^ /NDF 
of 5.0/6 and 13.3/6 at y = and y — 0.8, respectively. 
Near y=Q the antiprotons have a slightly smaller value 
of space-averaged phase-space density compared to that 
of protons, suggesting a small chemical potential. At 
V^i™ ~ ^'^■'^ GeV the antiproton space-averaged phase- 
space density was 38 times smaller [T^, which suggests 
a much larger baryo-chemical potential at the lower en- 
ergy. The inverse slope derived from the space-averaged 
phase-space density is T = 930 ±110 MeV for protons. 
This is consistent with data at ^s^^ — 17.3 GeV, but 
much higher than the value extracted at ^s^^ — 4.9 GeV 
where the inverse slope is about 350 MeV, (see Fig. [6]). It 
should be noted that the proton, kaon and pion spectra 



at midrapidity can be well described by blast-wave fits, 
which suggest that this increase in the inverse slope with 
Js^^ is largely driven by an increase in radial flow. This 
is supported by the fact that the phase-density of pions 
is characterized by a much smaller inverse slope (~ 140 
MeV) than that of protons [13] . 

Figures |4] and [5] imply that the volume of homogeni- 
ety, ^ l/i?2, and the space-averaged phase-space density 
vary little with rapidity at ^fs^ = 200 GeV. This is 
in stark contrast to the situation at lower energy. Ta- 
ble [n] shows the rapidity dependence of B2 and the 
space- averaged phase-space density at py = for cen- 
tral Pb+Pb collisions at ,/^ = 17.3 GeV [El [25l [27] . 
These data are compared to BRAHMS results at Pt/A 
— 1.3 GeV/c which show very little rapidity dependence. 
Finally, Fig. [6] shows the evolution of the space-averaged 
phase-space density near midrapidity as the energy of 
the system grows from AGS to RHIC values. The space- 
averaged phase-space density of protons decreases with 
energy while that of antiprotons increases. At AGS en- 
ergies (/) at rriT = is 10 times bigger than its value 
at RHIC, while at SPS it is 2 times bigger. In contrast, 
the pion average phase-space density has the opposite 
behavior as the energy of the colliding system increases; 
the values at RHIC [39l are higher than those at SPS by 
a factor ~ 2 where (/) = 0.45 [55j. At ,Js^A.^ GeV (/) 
was extracted from HBT and spectra measured at high 
rapidity (y) ~ 3.1 in Au+Au 10% central collisions and 
its estimated value at midrapidity is ~ 5 times smaller 
that the RHIC value [52]. For antiprotons the growth in 
(/) is similar to that of pions. The different trend for the 
protons may be driven mainly by baryon transport; at 
AGS, most of the beam protons end up at midrapidity, 
less so at SPS, and at RHIC, less than 10% of the beam 
protons are transported to midrapidity. The proton ra- 
pidity loss adds beam protons to the Gaussian rapidity 
density of produced protons and transforms the overall 
proton rapidity density into a flat distribution. 
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TABLE IL The rapidity dependence of B2 and the / for 
(Top) Pb±Pb collisions at ^/s^ = 17.3 GeV and (Bottom) 
^/s^ = 200 GeV AuAu coUisions. For the 17.3 GeV data 
the centrality is 23% at y=0.2 ^ and 20% at j/=0.8 and 1.3 

[HIESI. 



O 

o 



-a 

H 



10 

1 
10-' 
10-^ 
10-^ 
10- 

10-= 

10 

1 

10-' 

10-- 

10-' 

10--* 

10-= 




•—, 



••-... 



oo, 



■oo 



Oo. 



% 



+-f 



1.5<y<2.5 




0.5<y<1.5 



Oo, 



•Oo, 



'OOr 



^o^ 



°Or 



-f- °o o. 



2.8<y<3.2 



^ oO <? 



I ■ ... I ■ ... I ■ ... I 



0.5 1 1.5 2 2.5 3 3.5 

P^ [GeV/c] 



0.5 1 1.5 2 2.5 3 3.5 

P^ [GeV/c] 



FIG. 3. (color online) (anti-)proton and (anti-)deuteron pr spectra at various rapidities for the top 20% most central Au+Au 
collisions. The filled symbol part of the spectra show the pr intervals used in the coalescence analysis. Note the limits for the 
deuteron range are exactly twice those used for the protons. The errors are statistical only. The horizontal bars represent the 
bin width. 
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FIG. 4. (color online) B2 versus transverse momentum per 
nucleon at several rapidities for central Au+Au collisions at 
,Js^j^ = 200 GeV. Protons are displayed with filled circles and 
antiprotons with filled triangles. The solid line in the top left 
panel is an exponential fit to the data at y=0. This same line 
is shown, in dotted form, in the other 3 panels. The errors 
are statistical only. 



IV. SUMMARY 



The rapidity dependence of the deuteron production in 
— 20% central Au+Au collisions at 200 GeV has been 
studied in the context of coalescence models. Near cen- 
tral rapidity the proton and antiproton phase-space den- 
sities are very similar, suggesting a small baryon chemical 
potential. The coalescence parameters, or B2-values, for 
deuterons and antideuterons are also very close, suggest- 
ing similar freeze-out volumes for protons and antipro- 
tons. B2 increases with px as expected for a system 
undergoing transverse flow; flow introduces a correlation 
between position and momentum. At a given pT the 
deuteron coalescence parameter B2 is independent of ra- 
pidity, which would imply that the homogeneity volume 
for protons of a given pT is almost constant from y=0 to 
y=3. Baryon transport in these colliding systems may 
be affecting the rapidity density of protons. It is inter- 
esting to note that the radial flow varies weakly with 
rapidity |54j . The weak dependence of radial flow on ra- 
pidity offers an explanation of why the pr dependence 
of B2 does not depend on rapidity. It does not how- 
ever address the fact that the overall magnitude of B2 is 
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constant from y=0 to y=3. The proton space-averaged 
phase-space density also shows no significant rapidity de- 
pendence, whereas it depends strongly on ttit and y/s^. 
The space-averaged phase-space density can be thought 
of as the ratio of the number of protons per unity rapidity 
over the volume of the proton source. The rapidity den- 
sity of protons decreases by only a factor of 0.84 ± 0.01 
from y=0 to y=3. The constancy of the space-averaged 
phase-space density in rapidity is consistent with the sim- 
ilar behavior found for B2. Thus we observe a striking 
invariance of the proton source over several units of ra- 
pidity. 



FIG. 5. (color online) The (anti-)proton space-averaged 
phase-space density {f){mT/A) for central a/s^ = 200 GeV 
Au-|-Au collisions at several rapidities. The solid line in the 
top left panel is an exponential fit to the data at y=0. This 
same line is shown, in dotted form, in the other 3 panels. The 
errors are statistical only. 
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FIG. 6. (color online) The proton space-averaged phase-space 
density for central collisions as a function of ,Js^ and rriT 
at j/=0 [1311231 1^127) . The errors are statistical only. 
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